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Abstract 

We describe a general N-solitonic solution of the focusing NLSE in presence of condensate by the use of 
dressing method. We separate special designated class of "regular solitonic" solutions do not disturbing 
phases of the condensate at infinity. All regular solitonic solutions can be treated as localized perturbations 
of the condensate. We found an important class of "regular" solutions which is a small perturbation at 
certain moment of time. They describe the nonlinear stage of the modulational instability of condensate. 
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1. Introduction 

This research was motivated by intention to develop an analytic theory of freak (or rogue) waves in ocean 
and optic fibers. In the recent time the simplest and most universal model for description of these waves 
is the focusing NLSE (Nonlinear Schrodinger Equation). In application to the theory of ocean waves this 
equation is used since 1968 [l[. In nonlinear optics it was known even earlier Q|. 

The focusing NLSE is the model of the first approximation. For the surface of fluid this model describes 
the essentially weakly nonlinear quasimonochromatic wave trains with maximal steepness not more than 
0.15 only Q. In nonlinear optics its application is also limited to the case of waves of small amplitudes (see, 
for instance Q). Nowadays a lot of models generalizing the NLSE are developed. For the surface waves they 
are Dysthe equations for the waves in optic fibers equations include the third time derivatives and 

more complex forms of nonlinearity (see for instance 0, dj). Also, freak waves in the ocean were studied by 
numerical modeling of exact Euler equations for potential flow with free boundary 0, [l(| • The behavior of 
freak waves studied by NLSE and by more sophisticated models shows considerable quantitative difference. 
Nevertheless, advanced improvement of NLSE does not lead to any qualitatively new effects. That means 
that a careful and detailed study of NLSE solutions is still very important problem. 

There is a point of common agreement that extreme waves on deep water appear as a result of modu- 



lational instability of quasimonochromatic weaky nonlinear stationarity waves [9j, llll |12| . In terms of the 
NLSE model it means that we should study instability of the " condensate" in the focusing version of this 
equation. 

It is known since 1971 that the NLSE is a system completely integrable by the Inverse Scattering Method 
(ISM) [lj|. Since this time hundreds of papers and several books have been written to this subject (see 
for instance the monographs 14, lj| lj|. In this sea of literature one can find some articles devoted to 



development of the ISM for NLSE with focusing in presence of condensate (see, for instance. This ISM 
technique in application NLSE not is novelty 

However the central question appearing in this theory - what is long-time consequence of the modulational 
instability it is not answered so far. In this article we will study only evolution of a special class of localized 
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initial data presented by exact n-solitonic solutions of NLSE. Solitonic solution in presence of unstable 
condensate have a long history. 

In the nonsolitonic case a solution of auxiliary linear Zakharov-Shabat system (wave function) is analytic 
in the right half-plane of spectral parameter A = Xn + iXj with exeption of a cut on the real axis < \r < A, 
where A is the condensate amplitude. Each soliton adds a simple pole to some point of the right half plane 
including the cut. 

The simplest solution was found by E.Kuznetsov in 1977 [ljj. Later on this solution was rediscovered by 
other authors 1|| 18j]. In this case the cut is posed on real axis outside the cut Xji > A. The Kuznetsov 



solution is a localized bump oscillating in time. A period of oscillation growth as as \r — > A and become 
infinite in the limiting case \r = A. In this limit Kuznetsov's solution turns to the instanton found in 1983 
by Peregrine [20[. This is a pure homoclinic solution. It starts at t —oo from the condensate and return 
at t — > +oo to condensate with the same phase. The importance of these solutions for the development 
of the theory of freak waves is stressed in the article [25j . Later on multiinstanton solutions were found 

[HQ. 

In 1986 N. Akhmediev and V. Korneev discovered solitonic solution periodic in space but localized in 
time (Akhmediev's soliton) [27|. This solution is "almost homoclinic". It starts from a condensate and 
return to condensate with the same amplitude and different phase. In the Akhmediev case the pole is located 
inside the cut < Xjj < A. 

In a general case the pole is located in any point on the right half-plane. A generic solutions moves 
and oscillates. At x — > ±oo it goes to two condensates with the same amplitude and different phases. This 
solution was firstly found by M. Tajiri and Y. Watanabe in terms of Hirota's method. They also obtained 
general two-solitonic solution and suggested N-solitonic solution expression in 1998 [28]. However, in this 
article the properties of solitonic solutions except their stability were not studied. Later this solutions was 
obtained by some other methods and discussed in the frame of freak wave theory by A. Slunyaev at al. and 
N. Akhmediev at al. l3ll I33I. l34j|. In 2011 general one-solitonic solution was obtained by the authors of 



this article with the help of d problem [35j . 

In the recent time two-solitonic solutions and their degenerations were studied in the papers 28, In 
this article we describe a general N-solitonic solution of the focusing NLSE in presence of condensate. We do 
not insist on our priority in this point but we believe that our method for its construction is the most simple 
and economic. But this is not a major point of our article. This article is more "practically oriented". We 
give a partial answer to a major question - what is the nonlinear stage of the modulational instability. When 
we speak about this subject we must clearly separate development of two types of initial perturbations. One 
can consider periodic or quasiperiodic in space perturbations. Or one can study development of perturbation 
localized in space. The time-behavior of these two types of perturbations is different even in the linear theory. 
The nonlinearity amplifiers this difference enormously. In this article we will speak only on nonlinear behavior 
of spatially localized perturbations. We will show that a certain class of small perturbations can be described 
by 2N - solitonic solutions. Mare exactly speaking - by 2N "near Akhmediev" solitons. 

A general N-solitonic solution tends at x — > ±00 to the condensate with opposite phases. In this article we 
separate special designated class of "regular solitonic" solutions do not disturbing phases of the condensate 
at infinity. All regular solitonic solutions can be treated as localized perturbations of the condensate. In 
general case they never are "small". However we are able to construct an ample class of regular solitonic 
solutions which is a certain moment of time (suppose at t = 0) are small perturbation of condensate. In fact 
they are pairs of " near- Akhmediev" solutions. In this part of our article we study only the simplest case 
of 2-solitonic solutions belonging to this remarkable class. In the second part of the paper we discuss the 
general situation. 

2. NLSE via dressing method 

We study solutions of the following NLSE 

i V t-\^ x -M 2 -\A\ 2 )^ = i) (1) 
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with nonvanishing boundary conditions |</j| 2 — > \A\ 2 at x — > ±00. Without loss of generality we assume A is 
a real constant. Equation ([1} is the compatibility condition for the following overdetermined linear system 
for a matrix function "J [131 ]: 

i— = (\U + W)V (3) 



Here 



at 



7 = 



From ©, © one get 



9a; 



1 \ _ / (p 
-ir u_ v-y>* 



(5) 



(9/- 



= -*- 1 (AC7 + W) (6) 



and 



ax 

i^- = -*+(A*C/+ + iy + ) (7) 

We consider equations Q in the point A = —A*. One can see that 

U+(-X*) = -17(A) 

W + = W (8) 

Hence systems (0 , coincide. It means that they have a class of solutions "J satisfying the constrain. 

*+(-A*) = * _1 (A) (9) 

Thereafter we will assume that this condition Q is satisfied. 

The idea of the dressing method is following 36] . Suppose we know some solution 930 of the NLSE (fTJ) 
together with a fundamental solution V^o as a matrix function on x, t, A, satisfying condition of corresponding 
linear system. 

ox 

C/o and Wo are obtained by replacing ip — > tpo in Q, ([5]). Then we introduce "the dressing function". 

X = **o 1 (11) 
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(notice that \1/ is still unknown!). We demand that X is regular at infinity 



X (\)->E + Z + --- at |A|^oo 

*-(;;) (-) 

Apparently 

X + (-A*) = X- 1 (A) (13) 
The dressing function X satisfies to overdeterminated system of linear equations. The first equation is 

^ = U X ~xUo (14) 

and 

Ox' 1 
dx 

Equation (pT5|) can be rewritten as follow 



= - X - 1 U + U oX - 1 (15) 



U = -x(^ ~ Uoh" 1 (16) 

Now if choose dressing function x such that U defined from (I16p has no singularities on the A-plane we 
construct new solution of equation @ . 

According to the Louiville theorem in this case function U must completely be defined by its asymptotics 
at A — > oo. By plugging ((?]) to (|T^|) we find "dressing formula". 

u = u + [x,I] (17) 

or 

<P = fo - 2x ( i2) (18) 
So far we performed dressing equation @. From equation ((3]) we can derive following relation 



XU + W = -x (i^ - XU - Wo^x' 1 



(19) 



If we demand now that U, W has no singularities on all A - plane, including infinity (another word, we 
require that W does not depend on A), we realize "dressing" of equation §5§ either. After performing the 
dressing procedure and determination of X function 

* - A*o (20) 

satisfies equations ((2J, © where ip is given by (|18|) . This is a new solution of equation ((TJ. 

There several methods for constriction of dressing function \- The method based on the Rieman-Hilbert 
problem was described in article of Zakharov and Shabat published in 1979 36] More advanced approach is 
based on the use of the local 9-problem. It will a subject of our next paper. In this article we use "dressing 
for poor" making possible to construct multisolitonic solutions of NLSE on an arbitrary background. 
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3. General N-solitonic solution 

In this chapter we construct solutions of the NLSE following the method developed by V. Zakharov and 
A. Mikhailov [38|. Let us assume that ^o(x,t, A) is known and will consider that x is a rational function 
of A. 

x = ^+Err ( 21 ) 

- A-A„ 

Without loss of generality we can assume ReX n > 0. As x satisfy the condition (fTB")) . x -1 can be written as: 

n n 

Let us denote 

. .. U, 

Xn=X 



A=-A* 



^-EiT^V (23) 



and consider identity 

XX" 1 = E (24) 

in the point A = —A* identity (IM1) means that 

XnU+ = (25) 

In the same way 

xtU n = (26) 

Hence U n , U£ are degenerated matrixes, x„, Xn are degenerated either. Now we introduce two set of complex 
valued vectors p na , q na and put 

Una/3 = Pna<lnl3, ^nafi = InuVnp, 

Condition (|23|) means that 

X»<£a = (27) 

X«p = 5 a p + r- r 
» A _ A » 

x a0 = s«p-22 ttk (28) 

Let us plug (12"51) to equation (ITBl . In a general case function U acquires poles at points A = A„ and A = — A* . 
To perform dressing we must annulate residues in all poles. We consider equation (1161) in neighborhood of 
the point A = — A* . To kill the residue in this point we have to demand 

Xn - U (-X*S\ CaPlfS = (29) 

or 

xn(^-tM-A;)< Q + c^) =o (so) 

However in virtue of (f2"T| the last term in (j3"0)) is cancelled. Hence it is enough to demand 

^ - Uoi-KUq^ = (31) 
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This equality is resolved as follow 

Ina = *0,a/3(-A*)^„/3 

Here £ Q is arbitrary complex-valued vector. Thereafter we will use notation 

Fna/3 — *0,a^(— A*) 



(32) 



(33) 



Since this moment we consider that the set of vectors g„, ~<f * is known. To find the second p„, we 
must resolve relation (|27[) which is equivalent to the system of linear algebraic equations. 



m " L 

Here (q n ■ ~i* m ) — <Zn,i<Zm.i + Qn^Qm i is a scalar product of ~(f n and q m vectors. Let us denote: 



(34) 



M — 



A n + A* 



and M = det{M nm ). M is a Hermitian matrix: 



Now (j3"2")l can be rewritten as 

m 

We need to find \ from asymptotic expansion of \ (H21), which can be represented as: 

Xa/3 = ^2Pn, a qn,/3 

n 

This sum can be calculated as determinant ratio: 



(35) 



(36) 



M, 



Xa/3 



OL0 



M 



(37) 



Here M a p is extended matrix: 



M, 



a/3 



q\,p ■ ■ ■ q_n.fi 



QL 



Ml 



(38) 



Here M T - transposed matrix. We will find the the solution of NLSE ((T|) from condition (fT5|) as: 



M- 



4 12 



(39) 



Notice that formula (|3"T|) (f3"5 ]) in the simplest case dressing on zero background were found by Faddeev and 
Takhtajan [Hj]. 

Notice also that function \af3 also can be presented as ratio of two determinants. 



M, 



Xa/3 



a/3 

M 



(40) 



Here 



M afs = 



5 a 

Ql,a 
9n,a 



91,3 
A-Ai 



A — Ay, 



M 1 



(41) 



The dressing formula (|18[) can be written in the explicit form. Let 

^011 1*012 



*o = 



1 



021 1*022 



Then 



q, = JL ( ^11*011 + Mi 2 *o2i M n *oi2 + Mi 2 * 22 
X M [ M 2 i*oii + M 22 fern M 2 i*oi2 + M 2 2#o22 



Pn -> Pn 

a„. 



Now mention that transformation 

q n a n q n 

where a n are arbitrary complex numbers does not change the result of dressing. Thus one can put 

^ = I Cn 



(42) 



(43) 



(44) 



(45) 



Finally constructed n-solitonic solution depends on 2n sets of complex numbers A n , C n or of An real param- 
eters. Thereafter we will assume ReX > In fact, it is enough to enlist all possible solitonic solutions. This 
fact will be proven in a separate paper. So far we are sure that X n + A* re =/= 0, thus equations (|34p always are 
resolvable. 

Now we present one-solitonic solution. Function \ nas only one pole at A = 77, while x _1 has a pole at 
A = — rf . They can be presented in the following form (see [13j). 



X = E 



U 



X-r, 
U+ 
X + r]* 



As before 

U a /3 = Pat 

Vectors p a ,q a are connected now with relation 



Po 



As a result XiX 1 are 



Here 



X = E 



x- v ' 



{rj + r]*)P^ 
X + rj* 



a/3 



1a<lP 



ml 



1 92 1 



(46) 



(47) 



(48) 



(49) 



(50) 



P 2 = P thus P is projective matrix. As before we consider that the "seed" matrix ^o(x, t, A) is known. 
According to our defenition (j33|) 



Fa0 = ^0,ap(~V*) 
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(51) 



now 



<Zi — Fn + CF12 

q* 2 = F 2 i + CF 22 (52) 
C is arbitrary complex constant. Finally a new solution of NLSE is 

<P =l P0- -j 12 , I |2 (53) 

lair + M 2 

This formula presents the solitonic solution on an arbitrary background. It was established first in 1979 by 



Zakharov and Shabat [36] and reobtained in 1988 by Its, Rybin and Sail [39j (see also |41|). 

All results of this chapter can be extended to much more general class of nonlinear wave systems which 
can be presented as compatibility condition for overdeterminatcd linear system 

tf* = U(X)V 

* t = V(A)* (54) 
Where U, V are rational matrix 2x2 functions on A satisfying the follows condition 

U + (-X*) = -U(X) 

V + (-X*) = -V(X) (55) 

In particulary it is extended to all higher members of the NLSE hierarchy. One need only minor generaliza- 
tion to extend this procedure to the case of n x n matrix systems. 

4. AT-solitonic solution on condensate 

Since this moment we study dressing only on the condensate background. Now one can put tpo = A. 

U Q =( \ A x ) W Q = (56) 



And ^0 can be found as 



(x,t,X) . e -<j>(x,t,X) 

v/l-s 2 (A) V S ( A ) 



*oOM,A) = 7= == ( „, x x p -4,{x,t,x) ) ( 57 ) 



Here 



kx + flt, k 2 = X 2 -A 2 . fl = -iXk, s 



X + k 



Thereafter we assume that function k(X) = \/X 2 — A 2 has a cut —A < ReX < A. Thus fc(A) — > X at A — > 00. 
Then 

_ 1 If e -<t>{x,t,\) —s(X) ■ e~^ x,t ' x ^\ 

(X,t,X) = TT^lAT \-s(X) ■ e^'V e^'V ) (58) 

Notice that 

fc*(-A*) = -fc(A), s*(-X*) = -s(X), cj>*(-X*) = -4,(X) 

One can check that 

*o X ("A*) = *o + (A) (59) 



Thereafter we denote for simplicity. 



and by virtue of ((5 



4>n = 4>n(K) 



i{-K) = -<Pr, 



s(A„) 



(60) 



(61) 



Then 



*o(-a;) 



Q*n=F n 



1 

C n 



(62) 



(we can omit factor -. 1 because, as we mention before, it does not change the result of dressing)and 



Qni = e 



q n 2 



S n € 



Cle 4 



(63) 



So far we assumed that % is a rational function on the A-plane with cut —A < ReX < A. Now we perform 
the Jukowsky transform and map this plane onto the outer part of the circle of unit radius (see Fig. [1]). 



A = ^4) 



(64) 



The A-plane is mapped onto exterior of the circle of radius R. If the pole is posed at A = A„, A n = ■jiS.n + j-) 



Im(A) 



Im(£) 



-A 



Re(A) 




Re(£) 



Figure 1: Uniformization of the plane of spectral parameter with the help of Jukovsky transform. 

In new variables 



qm = e + — C*e 



q n 2 = j-e 



Finally we use parametrization 



Clef 



(65) 
(66) 



f — 7? p % 



(67) 
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and thereafter wc denote R n = e z ™. Now 



A 



Rn 

A 



iA 



A„ = — LR„ + — cos(a„) + — [R n - — sin(a n ) 



cosh(z„) cos(a„) + i sinh(z„) sin(a„) 



and we also denote 



Now 



w„ 



-ia„ 



e e 



to* 



Expressions for lf„ can be simplified by redefinition of phase factor <p n . One can put 



q n i = e ^ + e 



q n 2 = e 



Now 



Here 



Here 



u n + iv n 



^n cC n % *Ynt ~\~ ^ f-^n ^n — k n X LO n t <2^ n 



A 
2 
A 
~2 



se n = —[R n - — ) cos(a n ) = Asinh(z„) cos(a n ) 



Rn 



Rn 



In = -- 



sin(a„) = Acosh(z„) sin(a„) 
A 2 

R n + TT7 ) sin(2a n ) = — — cosh(22 n ) sin(2a„) 



Rl 



A 2 



A 2 



~[ R n~ ^2 ) cos(2a„) = — sinh(2z„) cos(2a„) 



The N - solitonic solution is invariant with respect to shifts in time and space. If one replace 

X — > X + x t — > t + t 

then 

fJ-n -> A*„ + 2(a3„a;o + 7„i ) + 2(fc„.x - w n *o) 

It means that a space-time shift leads to renormalization of constants C n . 

C \ C „2(aQ„2;o+7n*o)+2i(fc„a;o— v n to) 
y-^n * ^n^ 



(68) 
(69) 
(70) 

(71) 
(72) 

(73) 



(74) 

(75) 
(76) 
(77) 



7V-soliton solution can be considered as a nonlinear superposition of n separate solitons. Each of them is 
characterized by the group velocity. 



V Gr 



In 



(78) 
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and phase velocity 



V Pk~ = 'IT 



(79) 



In the case when all group velocities coincide the solution presents the "solitonic atom" first mentioned in 
[36| - a complicated configuration of solitons moving together. Complicated solitonic atoms containing a 
large amount of solitons should be described by method of statistical mechanics. In two-solitonic case 



-2 ^ sin(ai) 



R\ + 1 



• sin(a 2 ) 



Rx,i(Rl-V)~~'^' R2{Rl-l) 

Shapes of a typical "two-solitonic atom" and " three-solitonic atom" is presented on Fig. [5] and Fig. [3] 



(80) 



Re(<p), Im(p), M 2 ,A=1 




Re(<p), Im(^), \<p\ 2 ,A=1 




Figure 2: Two-solitonic atom at different moments of time. Ri = 1.5, R2 = 2,«i = ^,02 = 0.735242 Re(tp) - Green, Im(ip) 
Red, \ip\ 2 - Blue 



Re(^), Im(^), M 2 ,A=1 Re(^), Im(^), M 2 ,A=1 



25 r 25 




Figure 3: Three-solitonic atom at different moments of time. Ri = 1.5, R2 = 2,i?3 = 2.5«i = ^ ,«2 = 0.735242,03 = 0.6721 
Re(ip) - Green, Im(ip) - Red, \ip\ 2 - Blue 

If all group velocities are different the N - solitonic solution asymptotically at t — > ±00 separates to 
superposition of far remote from each other individual solitons. This fact makes possible to determine 
asymptotic properties of 7V-solitonic solution. We will discuss only the case when all z n > 0, thus the 
solution not contain "Akhmediev components". Now at x — ¥ ±00 

<p ->■ Ae la± (81) 
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Apparently phases a are constants in time. For far separated solitons 

a + = 2(ai H h a„) a~ = -2(ai H h a„) (82) 

This fact holds even in a general case when some group velocity could coincide. 

We developed the analytical code making possible to calculate iV-solitonic solutions by the use of the 
" Mathematica" software and checked relations directly for n = 2,3. 

If we are interested by ./V-solitonic solutions localized in finite domain of space and not perturbing far 
remote condensate we must put 

ai + --- + a n = (83) 

We call this solution "regular solitonic solution" and thereafter study only this case. If we assume that the 
modulation instability develops from localized perturbation, only regular solution can be used as model for 
its nonlinear behavior. 

Among one-solitonic solutions only Kuznetsov-Ma and Peregrine solution are regular. In two-solitonic 
case we can construct a broad class of regular solutions. 



5. One soliton solution on condensate 

One soliton solution on condensate background can be obtained by implementation of the results of the 
last part of §3. This solution is defined only by one complex eigenvalue Ai = r\. Let us denote 

= a = s{ v ) (84) 

remember that by virtue of ([55)1 

<K-rf ) = -0* a(-if) = -s* (85) 

After uniformization 

qi = e~* - C**se q 2 = -se-+ + C*e* (86) 
and the general one-solitonic solution is presented by expression. 

1 + |c| 2 ie 2 (^*) + C*e 2 * + Cjhe 2 ^ 

ip = A-A(^ + n— : 7TT* TXZS (87) 

V ; i + |q 2 e 2 W+**) +C*^±| F e 2 ^ + C 1 £±f F e 2 ^ V ' 

After our standard parametrization: 

£ = Re ia , C = e i8+ ^, R = e z 

we obtain (Thereafter we call R - amplitude and a - argument of a solution) 

(f) = u + iv 

u — asx — ft v = kx — ut 

A 1 A 2 1 

-(-R- -^)cos(a), 7 = - — (R 2 + — ) sin(2a) 

k=^(R+^)sm(a), w = ^(ii 2 --^)cos(2a) 

The general one-solitonic solution depends on four scalar parameters R, a, 9, \i. Two of them 8, fi are 
responsible for shift in time and in space. Thereafter we denote R = e z . Then 

A 2 

se = Asinh(z) cos(a), 7 = — — cosh(2z) sin(2a) 

A 2 

k = ylcosh(z) sin(a), cj = — sinh(2z) cos(2a) 
12 



If we put [i = 0, 9 = 0, the one-solitonic solution can be written as follows 



A 



cosh(,z) cosh(2it) + cos(a) cos(2u) 
cosh(z) cos(2a) cosh(2w) + cosh(2z) cos(a) cos(2i>) + i ( cosh(z) sin(2a) sinh(2u) + sinh(2z) cos(a) sin(2t>) 



To obtain the general solution one must replace t — > t — to, x — > x — xq where to 
This solution is localized in space if R ^ 1. In this case asymptotics of 



— 1M 

2 7 > 



X 



2 k ■ 



tp -> -Ae Zllx at 

-lia 



-00 

<p -> -Ae~ 2lQ at n +oo 
|(/?| 2 = A 2 at n ±oo 



We see that solution (|8"8"j) has identical asymptotics at a; — > ±oo only in the case a = 0, when the pole 
is posed on the real axis. 

The position of pole defines a type soliton. Different possible position are plotted in the Fig. 0] We 

Im(£) 




Re(£) 



Figure 4: 

"Species" of one-solitonic solution. ■- Kuznetsov-Ma, v - Akhmediev A - General, *-near Akhmediev, 

Peregrine 

describe below all these "species" of solitonic solution. 

Let us put a = 0, R > 1, z > 0. This is exactly Kuznetsov-Ma solution. If the solution is centered at 
x = and 9 = it is presented by following expression. Then 



.4 



tp : 



cosh(z) cosh(2w) + cos(2w) 

Here 

u = A sinh(z)a; 
This solution is periodic in time. Period of its oscillations 

47T 



cosh(z) cosh(2u) + cosh(2z) cos(2w) + i sinh(2z) sin(2t>) 



A 2 

v = — sinh(2z)i 



Asmh(2z) 
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(89) 
(90) 

(91) 



T ^ oo if z^O, R^l 
T ^ if z^oo, fi^oo 



(92) 



In the limit R — ¥ l,z — >■ Kuznetsov-Ma soliton turns to the standard standing soliton of NLSE on zero 
background. A typical behavior of Kuznetsov-Ma soliton is presented on Fig. [5] In the limit R — > 1 T — > oo 




Figure 5: Kuznctsov soliton at the moment of minimum(left) and maximum(right) of its amplitude. R = 2. Re(ip) - Green, 
Im(ip) - Red, |i^| 2 - Blue 



the Kuznetsov solutions turns to the Peregrine homoclinic soliton presented by expression 

1 - lit 



ip = A 1 -4 



' 1 + 4a; 2 + 4i 2 



(93) 



This solution was found by Peregrine in 1983 [20( and reopened by Its, Rybin and Sail in 1988 [39(. In 
their article it was called "exulton". This solution is actually "instanton" appearing from the condensate 
and disappearing. Indeed if — > A at t — > ±oo. More general " multiinstantonic" solutions are found in 
papers 23 , 24 1 . It is a popular idea to treat this solution as a model for freak waves 2|| 26| . Multisolitonic 
solutions describe a special scenario of modulation instability development, when growing perturbation in 
the long run return back to condensate. We will show in this article that this is a very special scenario, 
badly unstable with respect to small deformations. In a more general case waves developed from small 
perturbation do exist infinitely long time. The Peregrine soliton is plotted on Fig. [5] It was reported that 
this solution was observed experimentally in water wave tank [22| and optical fibers [2~lj . 

Let us put now R = 1, z = 0. Again we will put 8 = 0, fi = 0. We obtain famous " Akhmediev breather" 
periodic in space and localized in time. 



-A 



cos(2a) cosh(2it) + cos(a) cos(2u) + i sin(2a) sinh(2u) 



here 



1 



ii 



cosh(2it) + cos(a) cos(2w) 



A 2 sin(2a)i v = A sin(a)cc 



(94) 



Akhmediev solution is homoclinic in a "weak sense". It has following asymptotic behavior at t — > ±oo. 

tp -> -Ae~ 2la at t 
Lp -> -Ae 2ia at t 



-oo 
-oo 
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Figure 6: Peregrin breather at different moments of time. Re(ip) - Green, Im(tp) - Red, \ip\ 2 - Blue 



Hence \A + \ 2 = \A | 2 , but A + ^ A .In the especially interesting case a = j A + = —A . Spatial period of 
the Akhmcdicv breather 



L = 



2tt 



A sin(a) 



(95) 



Tends to infinity if a — > 0. In this limit solution tends to a periodic set of distantly remote Peregrine solitons. 
The Akhmediev breather is plotted on Fig. [JJ 

Now we return to the general one-solitonic solution. In the explicit form it was found in the paper of 
Tajiri and Watanabe (2|| and reobtained then by several authors [2^, [33|, HO, [35| . This is a localized object 



of size / 



Acos(a)sinh(z) propagating along the condensate with group velocity 

cosh(2z) 



V„ 



1 

a; 



A- 



sinh(;z) 



sin(a) 



(96) 



Soliton's amplitude oscillates with angular frequency lo. Soliton is filled with carrying wave propagating 
with phase velocity. 

„..,,, cos(2a) , . 

V ph = 2Asinh (z) . ) / 97 

If a —> this carrying wave vanishes. A typical "general" soliton is plotted on Fig. [5J 

It is important to average by integration over the period of oscillations T = — the squared absolute value 
of onesolitonic solution (155)) . It demonstrates the behavior of average value of the Hamiltonian. Remarkably 
that the answer can be obtained analytically in very simple form. This requires move to a coordinate system 
which travel with group velocity V gr 



x + —t 

as 
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(98) 
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10 



Re((p), Im(<s), M \ A=l 

t II " 



-10 



-2 




10 



Re(^), Im^), \<p\ \ A=l 
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10 



Re(<p), Im(<p), M 2 ,A=1 

or 
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Figure 7: Akhmediev breather at different moments of time. H = l,a=-?,/Lt = 0. Re{ip) - Green, Im(ip) - Red, </}| 2 - Blue 



Then the dependence on time of hyperbolic functions disappear and they can be regarded as constant when 
integrating over time. As a result 



< M >t-- 



2tt 



2 /■27T 



1 



cosh(z) cosh(2it) + cos(a) cos(r) 



cosh(z) cos(2a) cosh(2u) + cosh(2z) cos(a) cos(r) 



+ 1 cosh(z) sin(2a) sinh(2w) + sinh(22:) cos(a) sin(r) 



eh 



(99) 



Integral is taken by residues after the standard change of variables w = e lT . The path of integration now 
| to | = 1. The integrand has one pole of the first order at point zq — and two poles of the second order in 



points Z12 = — 



cosh(2) cosh(2u) _j_ /cosh 2 (z) cosh 2 (2-u) 



cos(a) 



cos 2 fa) 



— 1. However, only wq and w± lie inside the circle \w\ = 1. 



After integration we need back to initial coordinate system. The final answer is 

4cosh(2w) sinh 2 (z) cos 2 (a)(sinh 2 (z) + sin 2 (a)) 



< \<p\ 2 > T = A [1 + 



1 3/2 



cosh (z) 



(100) 



V cosh 2 (2u) - 

v ' cosh 2 (z) 

The case from Fig. [5] averaged by time is presented on Fig. [5] It is interesting to perform averaging of 
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Re(<p), Im(<p), M 2 ,A=1 




Re((p), Im(y>), M 2 ,A=1 




-10 



Figure 8: General solution at the moment of minimum(left) and maximum(right) of its amplitude. R = 2, a = Re(tp) 
Green, Im(ip) - Red, |</j| 2 - Blue 

<M 2 >„A=1 

6r 



-5 



-2 L 

Figure 9: i? = 2, a = ^ 



Kuznctsov-Ma solution in the point x — 0. In this case a — 0, u = and 

l2 2 / sinh (2) 

<M >r=A ( 1+ ^h%y 

In particulary in the limiting case we obtain the expected result 

< M 2 > T = A 2 

Now we pay attention to the special case of general solitonic solution when R 
" near- Akhmediev" solution of very large size 

1 



2zAcos(a) 



moving with very high group velocity 



Vc 



A 



qroup • / \ 

sm(a)z 
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(101) 
(102) 

1, z -} 0. This is the 
(103) 

(104) 



and very low phase velocity 

2Azcos{2a) 

Vph « r-^ — (105) 

sm(a) 

The soliton has inner quasiperiodic structure with characteristic scale I ~ A ■ If one observe this 
solution from a fixed point (suppose being at x = 0) the total time of this soliton passing is 

- 1 ' (106) 



F sroup 2A 2 cos(a) sin(a) 



T — >• 00 if a — > 0. It is interesting that this time does not depend on z. It is important to stress that the 
" near-Akhmediev" soliton remain after its passing slowly decaying " tails" . The " near- Akhmediev" solution 
is plotted on Fig. [TO] 



lllllllllj! 1 !!)" 



-200 -100 100 200 

Figure 10: The "near-Akhmediev" solution. (R = 1.02, a = -S) 



6. Regular two solitonic solution 

A general two-solitonic solution on an arbitrary background can be obtained after application of the 
dressing procedure described in §2 in the relatively simple from. 

N 

V=A-2- (107) 



_ B\ B 2 B3 i B 



A = 



A2 + A* + A2 A2 + Ai Ai + A^ 

kiW (titSXtlfa) 



(Ai + A^)(A 2 + A5) (AJ + A 2 )(A* + A X ) 
Bi = \q 2 \ 2 qnqi2 B 2 = (glfftOgli Q12 B 3 = (tfi^D^iife #4 = |?i| 2 <?2i522 (108) 

Let us denote also function S 

S = gn<722 - 912921 (109) 
18 



One can check the two following identities are valid. 

B x - B 2 - B 3 + B A = 

|gi| 2 |g2| 2 - {qiqViiqiqV) = 1 911*22 - qi2q2i\ 2 = \S\ 2 

Thereafter we will study the regular two-solitonic solution such that 

a2 = — cti = —a 

Now 

cosh(zi) cos(a) + isinh(zi) sin(a) 



(110) 
(111) 

(112) 



Ai ' 7^) l ' l,M " : ' ■ y( Ri If ^ ' A 

A ( l \ , x iA( 1 , 
A 2 = — ( R 2 + — ) cos(a) — ( R 2 — — I sm(a) = A 



R-. 



R<. 



cosh(z2) cos(a) — zsinh(z2) sin(a) 



(113) 



To simplify the final expression for two-solitonic solution we perform corresponding shifts in time and space 
making \C\ \ = \C 2 \ = 1. Hence now 

d = e Wl C 2 = e l ° 2 (114) 

Finally the solutions depends on five parameters Ri,R 2 , a, &i,0 2 . Let us introduce combination of two last 
parameters. 

e+ = e l + e 2 o- = e 1 -e 2 (115) 

Now 



qn = e -0i + e -ia-Z!+4>i 
q 12 = g-ia-'i-'h + e <h 



q 21 = e-^ + e ia - Z2+ ^ 
q 22 = e ia - Z2 -* 2 + 



(116) 



(pi = u\ + ivi 

U\ = 3d\X — -fit 

u 2 = ac 2 .x - j 2 t 



02 = u 2 + iv 2 
vi = kix — uj\t — 

v 2 = k 2 x — u> 2 t - -62 

aei = — f Ri — ] cos(a) = ylsinh(zi) cos(a) 

2 V R iJ 

ki = — -I Ri + — I sin(a) = ylcosh(zi) sin(a) 
2 V Ri J 

A 2 ( 1 \ A 2 
71 = -— f i?? + _ J sin(2a) = — — cosh^) sin(2a) 

A 2 f l \ A 2 
loi = —— I R\ — -2 ) cos(2a) = — sinh(2zi) cos(2a) 

4 V R \J 2 

ae 2 = — -( R 2 — —J cos(a) = A sinh(z 2 ) cos(a) 
2 V R 2 J 

k 2 = ~^{r2 + jj-j sin(a) = -Acosh(z 2 ) sin(a) 

A 2 ( 1 \ A 2 

72 = — f Rl + ^2 j sin(2a) = — cosh(2z 2 ) sin(2a) 

A 2 ( 1 \ A 2 
lu 2 = — ( R 2 - -^2 J cos(2o;) = — sinh(2z 2 ) cos(2a) 

19 



(117) 



In the special symmetric case 



A2 — Aj 
as 2 = aei lo 2 = ui 
k 2 = -ki 72 = -71 (118) 



and solution can be written in an direct analytic form 



here 



A M + iK , s 

tp = A-2A — (119) 

H 



H = 4 [cos 2 (a) + sinh 2 (z)] [cosh 2 (z) sin 2 (a) cosh(4a8x) — sinh 2 (z) cos 2 (a) cos(4fcir, — - )] 

+ sinh 2 (2z) cosh(47t) - sin 2 (2a) cos(4wt + 0+) 
+2 sinh(2z) sin(2a) sinh(z) sin(a) x 
[cosh(2ascc - 2jt) cos(2fcz; + 2wt + 2 ) + cosh(2aea: + 2jt) cos(2fcr - 2ut - 0i)] 

+2 sinh(2z) sin(2a) cosh(2:) cos(a) x 
[sinh(2fBx - 2-/t) sin(2fcx + 2cot + 2 ) + sinh(2fBa; + 27*) sin(2fcx - 2cot - 6{)] (120) 



M = sinh(2z) sin(2a) x 
^sinh(2z) sin(2a) [cosh(47t) + cos(4wi + 0+)] 

+2 sinh(z) sin(a) [cos 2 (a) + cosh 2 (z)] x 
[cosh(2aea; - 2jt) cos(2fcr + 2ojt + 2 ) + cosh(2eex + 2jt) cos(2fez; - 2ut - X )] 

+2cosh(z) cos(a) [sin 2 (a) — sinh 2 (.z)] x 

[sinh(2a3a; - 2jt) sin(2fca; + 2ut + 2 ) + sinh(2aex + 2jt) sin(2fc:r - 2ujt - X )] J (121) 



K = sinh(2z) sin(2a) x 
sinh(2z) sin(2a) sinh(47i) - sin(2a) sin(4a;t + + ) 

— 2 cosh(z) sin (a) [cos 2 (a) + sinh 2 (z)] x 
[cosh(2a8x - 27*) sin(2fcx + 2ut + 2 ) - cosh(2aex + 2jt) sm(2kx - 2ujt - 0i)] 

— 2sinh(z) cos(a) [cos 2 (a) + sinh 2 (2:)] x 
[sinh(2se:c - 2~ft) cos(2fcx + 2ut + 2 ) - sinh(2a8:r + 27*) cos(2fca; - 2ut - 0i)] (122) 

The regular two-solitonic solution in a general non-symmetric case and in the special symmetric case are 
plotted on Fig. [TT] and Fig. [T2] 
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Re(»>), Im(y>), \ip\ 2 , A=l Re(»>), lm(<p), \ip\ z , A=l 




Figure 11: Twosolitonic solutions. R = 2, a = 




Figure 12: Twosolitonic solutions. i?i = 1.5, ai = ^,i?2 = 2.5, 02 = —\ 



7. Annulation of solitons 

Let us consider the limiting case 

Now 



Rl — i?2 — 1 



Ai = = A2 = A2 = j4cos(a) 
aei = ac2 — wi = LJ2 = 
A 2 

k\ = -k 2 = Asin(a) 71 = -72 = — sin(2a) 
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Looking at ( UOSp we see that now 



N = — 1 r (B 1 - B 2 - B 3 + B A ) = (125) 
2A cos(a) 



til = — u 2 = — sin (2a)t 
Vi = Asin(a)x — u 2 = — Asin(a)x — ^#2 (126) 



Thus 



Here 



Now 



1 1 

-01 02 = —0 — -#2 



(127) 



A 2 

— sin(2a)t + iA sin(a)x (128) 



q n = e" 0+ 5 ei + e- ia +^-i e i g 21 = e^+i 02 + e™^^ 62 

gi2 = e-^-^+i 01 + e't'-i 91 q 22 = e la+4>+ i 02 + e ~<l>-i e ' (129) 

Vectors ~q*i,~q*2 & r e periodic functions on x and exponential functions on time. However simple calculations 
shows that A does not depend on x, t. Now 

/e+\ 

S = <?ii<722 - <7i2<72i = -4 sin(a) sin I — J (130) 



Hence finally 



sin 2 (a) . 2 / 9 + ', 



A 2 cos 2 (a) V 2 



The denominator A in this case does not depend on x and t. This is just a number. A 7^ if 9 + 7^ (we do 
not consider the special "Peregrine" case a = 0). Since iV = this means that in a general case 9 + ^ two 
pure Akhmediev solitons with opposite values of a completely annulate each other. In this case the dressing 
function \ is presented by the scalar matrix. 



A- Acos(a) J V 1 

In the case 9 + = 0, both numerator N and denominator A in (|131[) are zero. This case must be considerate 
separately. 



8. Small perturbation of condensate 

This chapter is the central in our paper. The remarkable fact of exact cancelation of two Akhmediev 
solitons with opposite values of angular parameter a makes possible to construct special and very important 
class of 2A^-solitonic solutions describing evolution of infinitesimally small perturbation of condensate. All 
of these solutions are regular. We start with the case N = 1 and assume that i?i, R 2 are close to one (see 
Fig. US). 



Ri ~ 1 + e R 2 ~ 1 + at 
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(133) 





1 £ 

— R=l+e 

/or 1 










\ / 1 

— R=l+ae 

^2 





Re(£) 



Figure 13: Poles location corresponds to small perturbation of condensate 



e > 0, a > 0, |e| <C 1 Deviation (5c/? of a solution from the condensate at the moment of time t — is given by 
expression 

~6N 



5(p 



(134) 



Here SN is a variation of numerator while A is a "modified" versions of the denominator calculated with 
higher accuracy. To calculate SN we noticed that 



Then 



Then 



Ai = — ( Ri H ) cos(a) + - — ( R\ I sin(a) ~ A cos(a) + ieA sin(a) 

2 V RiJ 2 V Ri. 

A 2 = \ Ri + -p— Jcos(a) — ^ ( i?2 — ^- )sin(a) « Acos(a) — iaeAsin(a) 
2 V R%J 2 \ i? 2 , 



Ai + A 2 = 2Acos(a) + + a)eAsin(a) 
A* + A 2 = 2 A cos(a) — i(l + a)eA sin(a) 



SN 



ill + a)e , . ,~ ~ . 
-i f sin a 5 3 -5 2 



AA cos 2 a 



(135) 
(136) 

(137) 
(138) 



Here Bi,B% are "modified" versions of B2,B^. To calculate A we assume that 9 + ^ and neglect by 
deviation of Ai, A2 from its limiting value Acos(a). Hence as before 



A = 



1 



AA 2 cos 2 (a) 



1911922 - 912921 1 



\S\< 



AA 2 cos 2 (a) 



(139) 



B 2 and B 3 are presented by expressions (|108p . Vectors ~<f 1, are given by formulae (|116l) . but to calculate 
"modified" functions we have to keep the first order of small parameter e in functions 4>\,4>2- After that 
functions <pi , <p2 should be considered as not depended on e in the small parameter expansion. Now 



k\ = —k2 = Asin(a) 
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(140) 



and 



Now 



Here 



aei = Ae cos(a) ae2 = aAe cos(a) 



h=4> + ef--6 1 



aef - -6 2 



<fi = iAsin(a)x 
Vectors q\, qi has following components 



/ = A cos(a)a; 



qil = e -4>-e/+§0l + e -ia+4>+tf~l;0l 
qi2 = e -i*-<t>-ef+i6i + e 4>+tf-¥i 



q 21 — e 4>-aef+±6 2 _|_ e ia-4>+atf- 
q 22 = gia+0-ae/+f 8 2 _|_ g -0+a£/- 



Lct us denote: 



Now 



* = Asin(a) ( = Acos(a) 



S = —Ai sin(a) 

Hence 

~ 4 sin 2 (a) 

A = — 

A 2 cos 2 (a) 

Then we noticed that 



V 2 



6+ 



cos — I sinh(e(l + a)(x) — i sin — I cosh(e(l + a)(x) 



[cos 2 — sinh 2 (e(l + a)(x) + sin 2 — cosh 2 (e(l + a)(x)] 



6+ 



B 3 - B 2 = (9u921 + 9i 2 922)<?2l'?12 - (911921 + 9l2922)9ll922 

= 922921 ( 1 9n I 2 - I912I 2 ) - 912911G921I 2 - I922I 2 ) 



Hence 



Sip = —2eAi(l + a) sin(a) 



B3 — B2 



(141) 
(142) 
(143) 

(144) 
(145) 

(146) 
(147) 

(148) 
(149) 



Function S are slow oscillating in space. The products (712922 and 9 2 i9*i consists both "fast" and "slow" 
components. Hence ip can be present as a sum of two terms 



Here 



5<pi — eAi(l + a)- 



5ip = Sipi + 5(f2 

sin(0+) 



cos 2 (^) sinh 2 (e(l + a)(x) + sin 2 (^) cosh 2 (e(l + a)(x) 
is " slow varying in spase function" , while 

1 

Sif2 = -eA(l + a) 



(150) 



(151) 



cos 2 (^) sinh 2 (e(l + a)(x) + sin 2 (^) cosh 2 (e(l + a)(x) 
sin(a) [sinh(2eaCa;) sin(2*a; - 6>i) + sinh(2eCx) sin(2*x + 2 )] 

+i cos(a) [cosh(2ea(» sin (2*2; - 61) - cosh(2e(>) sin (2*2; + 8 2 )] 
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(152) 



is fast oscillating. In symmetric case a = 1 and Sip becomes 



„ sin(a) sinh(2eCa;) sin(2*x - V) cos (^r) _ * cos(a) cosh(2eCx) cos(2^a; - V) sin (^r) 

Sip = —AeA /a . - ' „,., , 7. — (153) 

cos 2 (^-) sinh 2 (2eCx) + sin 2 (^) cosh 2 (2eCx) 



The most simple case appear when 6\ = 02 = \ an d a = 1 



Y cosh(2eCa;) v ; 

The typical small localized perturbation of the condensate at the moment t — is presented on Fig. [TU 
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Figure 14: The typical small localized perturbation of the condensate at the moment t = Da)o = l,e = 0.05, a = 
f , 6» 2 = f b) a = l,e = 0.05, a = f^, 0i = tt, 2 = c) a = 1, e = 0.05, a = ^, 6>i = vr, 6» 2 = d) a = 2, e = 0.05, a = 

7T /) 7T 



The comparison of exact solutions and its first order expansion (|f 34[) is presented on Fig. [15] 
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Figure 15: Exact solution (continuous curves), first order expansion (dashed curves) a = 1, t = 0.05, a = jjj, 9i = ^, 02 = \ 



We end up with the following result. The special class of centered in space and time two-solitonic 
solutions presenting. At the moment of time t = a small perturbation of condensate is described by five 
parameters a, e, a, 9 + ,0~ . Moreover this solution can be shifted in space and in time, thus the total number 
of free parameters is seven. In construction of such solutions one must be sure that 6 + ^ 0. The degenerated 
case 8 + — > must be studied separately. 

The derived 2-solitonic solution describe following physical process. 

Initially small localized perturbation of the condensate generates a pair of near-Akhmediev solitons 
propagating with very fast group velocity in opposite directions. In the symmetric case a — 1 these solitons 
are symmetric. Typical development of these small localized perturbation of the condensate is presented 
on Fig. [16J Fig. [TT1 Fig. [TH Fig. [19] In the second part of our paper we will show that this results can 
be essentially generalized. We construct 2N solitonic solutions presenting at t = small perturbation of 
condensate. Also in the second part of paper we will study specially degenerated cases 8 + = 
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